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ABSTRACT 

If double neutron star mergers leave behind a massive magnetar rather than 
a black hole, a bright early afterglow can follow the gravitational wave burst 
(GWB) even if there is no short gamma-ray burst (SGRB) - GWB asso ciation 



or th ere is an association but the SGRB does not beam towards earth (jZhang 



20121 ). Besides directly dissipating the proto- magnetar wind, we here suggest 
that the magnetar wind could push the ejecta launched during the merger pro- 
cess, and under certain conditions, would reach a relativistic speed. Such a 
magnetar-powered ejecta, when interacting with the ambient medium, would de- 
velop a bright broad-band afterglow due to synchrotron radiation. We study this 
physical scenario in detail, and present the predicted X-ray, optical and radio 
light curves for a range of magnetar and ejecta parameters. We show that the 
X-ray and optical lightcurves usually peak around the magnetar spindown time 
scale (~ 10 3 — 10 5 s), reaching brightness readily detectable by wide-field X-ray 
and optical telescopes, and remain detectable for an extended period. The radio 
afterglow peaks later, but is much brighter than the case without a magnetar 
energy injection. Therefore, such bright broad-band afterglows, if detected and 
combined with GWBs in the future, would be a probe of massive millisecond 
magnetars and stiff equation-of-state for nuclear matter. 
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Introduction 



The next generation gravi tational-wave detecto rs, such as Advanced LIGO (lAbbott et al. 



20091 ) and Advanced VIRGO (lAcernese et al.ll2008l ). are expected to detect GW signals from 
mergers of two compact objects. These gravitational wave bursts (GWBs) have well defined 
"chirp" signal, which can be unambiguously identified. Once detected, the GW signals would 
open a brand new channel for us to study the universe, especially the physics in the strong 
field regime. Due to the faint nature of GWs, an associated electromagnetic (EM) emission 
signal in coincidence with a GWB in both trigger time and direction would increase the 
signal-to-noise ratio of the GW signal, and therefore would be essential for its identification. 

One of the top candidates of GWBs is merger o f two neutron stars (i.e. NS-NS 
mergers) (ITaylor fc Weisbergl Il982t Kramer et al.l |2006[ ) . The EM signals associated with 



such an event include a short gamma-ray burst (SGRB) ( Eichlc r et al.lll989l ; iRosswog et al. 



2012 



Gehrels et al.l 



(ILi fc Paczvriski 



2005 



1998 



glow ( iNakar fc Piran 



Barthelmv et al 



Ku 



2011 



karni 



2005 



2005 



Berger 



Metzger et al. 



Metzger fc Bergerll2012l 



201lh . an optical "macronova' : 



20101). and a lon g lasting radio after- 



Piran et al.l 120121 ). Numerical simula- 



tions show that binary neutron star mergers could eject a fraction of the materials, forming a 
mildly anisotropic ou tflow with a typical velocity about 0.15 — 0.25c, and a typical mass about 



10" 4 ~ 1O~ 2 M fe.g. lRezzolla et al.ll201ll ; IRosswog et al.ll2012l ; iHotokezaka et al.ll2012f ). The 
radioactivity of this ejecta powers the macronova and the interaction between the ejecta and 
the ambient medium is the source of radio afterglow. Usually, the merger product is assumed 
to be a black hole or a temporal hyper -massive neutron star which suryives 10-100 ms be- 
fore collapsing into the black hole fe.g. IRosswog et al.ll2003t lAloy et al.ll2005l ; IShibata et al. 



20051 ; iRezzolla et al.ll201ll ; IRosswog et al.ll2012l ). Nonetheless, recent observations of Galactic 
neutron stars and NS-NS binaries suggest that the maximum NS mass can be h igh, which 
is close to the total m ass of the NS-NS system s (IDai et al.l 120061 ; IZhangl 12012k and refer- 
ences therein). Indeed, IHotokezaka et al.l (120121 ) show that for a wide range of parameters, a 
hyper-massive neutron star can be formed and survive for an extended period of time. Even 
for a very stiff equation of state for nuclear matter, as suggested by Dai et al. (2006) and 
Zhang (2012), a stable massive neutron star appears after the merger. This newborn mas- 
sive neutron star would be differentially rotating. The dynamo mechanism may operate and 
generate an ultra-stron g magnetic field (IDuncan fc Thompson! Il992j ; iKluzniak fc Ruderman 
19981 ; IDai &: Lulll998bl ). so that the product is very likely a millisecond ma gnet ar. Evidence 



of a magnet ar following some SG RBs has been collected in the Swift data (IRowlinson et al. 



2010l ; IRowlinson fc O'Brien! |2012| ). and magnetic activities of such a poster- merger massive 
neutron star have been suggested to interpret several X- ray afterglow properties of some 
SGRBs (foal et al.ll20o3 : iGao fc Fkrfood : En fc xJbooj ). 
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Recently, IZhangl (120121 ) suggested that if the merger product is a millisecond, massive 
magnetar, the magnetar wind can have a much wider solid angle than an associated SGRB, 
if any. Therefore, a bright X-ray afterglow due to dissipation of the magnetar wind can be 
associated with a GWB even if there is no SGRB-GWB association, or the SGRB does not 
beam towards earth. Numerical simulations suggest that a magnetar wind can be collimated 
by the ejecta to an opening angle of 30° — 40°. The magnetar wind can leak from such a 
funnel and dissipate through internal-collision induced magnetic reconnection and turbulence 
(IZhang fc Yanll201ll ). At a larger angle, the magnetar wind is confined by the heavy ejecta 
launched during the merger phase. The magnetar w ind can be d issipated upon hitting the 
ejecta, which again powers a bright X-ray afterglow ( )Zhangll2012l ). On the other hand, after 
this phase, a significant fraction of magnetar spindown energy would be used to push the 
ejecta, giving a stronger interaction with the ambient medium than otherwise (See FigJTJfor 
a cartoon picture for various emission components in this scenario). We study the dynamics 
of such an interaction in detail, and calculate the broadband afterglow emission from such a 
scenario. 



The model 



The postmerger hyper- massive neutron star may be near the break up limit, so that 
the total spin energy E rot = (1/2) 10% ~ 2 x IO^I^Pq^ er g (with / 45 ~ 1.5 for a massive 
neutron star) may be universal. Here Pq ~ 1 ms is the initial spin period of the proto- 
magnetar. Throughout the paper, the convention Q = 10 n Q n is used for cgs units, except 
for the ejecta mass M e j, which is in unit of solar mass M . Given nearly the same total 
energy, the spin-down luminosity and the characteristic spin down time scale crit ically depend 
on the polar-cap dipole magnetic field strength B p (IZhang fc Meszarosl 120011 ). i.e. L sd = 

10 49 erg s" 1 £2 DfPr 4 



£ S d,o/(l +t/T sd ) 2 , where L sdj0 



T sd ~ 2 X 10* S 745^15^0,-3^6 
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erg s " tfp j^itgi-g _2 3 , and the spin down time scale 



E ro t/ L : 



sd,0- 



Afte r the intern al dissipation of the magnetar wind that powers the early X-ray af- 
terglow JzhaneJ l20lA . the remaining spin-down energy would be added to the blastwave. 
The dynamics of the blastwave depends on the magnetization parameter a of the magnetar 
wind after the internal dissipation. Since for the confined wind, magnetic dissipation occurs 
upon interaction between the wind and the ejecta, in this paper, we assume that the wind 
is still magnetiz ed (moderately high a), so that there is no s trong reverse shock into the 
magnetar wind ( Zhang fc Kobayashi 2005 ; Mimica et al. 2009 1^1. As a result, the remaining 



If, on the other hand, the wind is already leptonic matter dominated, a reverse shock can be developed, 
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Fig. 1. — A cartoon picture of our model. A hyper- massive neutron star is formed at the 
central engine. Near the spin axis, there might be a SGRB jet. An observer towards this jet 
(red observer) would see a SGRB. At larger angles (yellow observer), a fre e magnetar wind 
may be released, whose dissipation would power a bright X-ray afterglow ( jZhangi 120121 ). At 
even larger angles (orange observer), the magnetar wind is confined by the ejecta (green 
shell). The interaction between the magnetar wind and ejecta would trigger magnetic dissi- 
pation of the wind and power a bright X-ray afterglow (jZhangil2012l ). After releasing some 
dissipated energy, a good fraction of the spindown energy would push the ejecta and shock 
into the ambient medium. Synchrotron emission from the shocked medium (red shell) would 
power bright X-ray, optical and radio afterglow emission, which is calculated in this paper. 
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spin energy is continuously injected into the blastwave with a luminosity Lq = £X s d,cb where 
£ < 1 denotes the fraction of the spin down luminosity that is added to the blastwave. The 
evolution of the blastwave can be describe d by a system with continuous energy injection 



f lDai fc Lulll998al : IZhang fc Meszarod 120011 ). 



The newly formed massive magnetar is initially hot. A Poynting fl ux dominated outflow 



is launched ~ 10 s later, when the neutrino-driven wind is clean enough ( jMetzger et al.ll201ll ). 
At this time, the front of the ejecta traveled a distance ~ 6 x 10 10 cm (for v ~ 0.2c), with a 
width A ~ 10 7 cm. The ultra-relativistic magnetar wind takes ~ 2 s to catch up the ejecta, 
and drives a forward shock into the ejecta. Balancing the magnetic pressure and forward 
shock pressure in the ejecta, one gets rest-frame sound speed v s ~ 10 _4 cLQ^ 2 7 A^ 2 M e j j^. So 

]/2 1/2 1/2 

the forward shock would cross the ejecta in t& ~ A/v s ~ 3 s L Q J 7 A/ M e /_ 3 . 

After the forward shock crosses the ejecta, the forward shock ploughs into the ambient 
medium. The dynamics of the blastwave during this stage is defined by energy conservation^ 

L t = ( 7 - l)M ej c 2 + ( 7 2 - l)M sw c 2 , (1) 

where M sw = ^-R 3 nrrip is the swept mass from the interstellar medium. Initially, (7 — 
l)M e jc 2 (7 2 — l)M sw c 2 , 7 — 1 of the ejecta would increase linearly with time until t = 
min(T sd , T dec ), where the deceleration timescale T dec is defined by the condition (7 — l)M e jC 2 = 
( 7 2 — l)M sw c 2 . By setting T dec ~ T sd , we can derive a critical ejecta mass 

M ej)Cjl ~ l^M & n^I%%fp-% 2 e /A ~ ^M^ltB^R^P^ 7 '*, (2) 

which separate regimes with different blastwave dynamics. For a millisecond massive mag- 
netar, the parameters J45, Rq, Po,-3 are an essentially fixed values. The dependence on n 
is very weak (1/8 power), so the key parameters that determine the blastwave parameters 
are the ejecta mass M e j and the magnetar injection luminosity Lq (or the magnetic field 
strength B p ). If M e j < M e j c l (or T dec < T sd ), the ejecta can be accelerated linearly until the 
deceleration radius, after which the blastwave decelerates, but still with continuous energy 
injection until T sd . Conversely, in the opposite regime (M e j > M e j jCj i or T sd < T dec ), the 
blastwave is only accelerated to T sd , after which it coasts before decelerating at T dec . In the 
intermediate regime of M e j ~ -^ej,c,i (or T dec ~ T sd ), the blastwave shows a decay after being 
linearly accelerated. 



which would predict additional radiation signatures (|Dail 120041 ) . 



2 The accurate expression for EqU] should be L t = (7 — 7 c j !0 )M c jC 2 + (j 2 — l)M sw c 2 , where 7 c j,o is the 
initial Lorenz factor for the ejecta, which we take as unity for convenience. 
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There is another critical ejecta mass, which defines whether the blastwave can reach a 
relativistic speed. This is defined by -E rot £ = 2(7 — l)M e jc 2 . Defining a relativistic ejecta as 
7 — 1 > 1, this second critical ejecta mass is 

M ej , c , 2 ~ 6 x 10~ 3 M Q h 5 P %C (3) 

An ejecta heavier than this would not be accelerated to a relativistic speed. 

Below we discuss four dynamical regimes. 

Case I: M e j < M e j jCj i or T sd > T dec . This requires both a small Lq (or low B p ) and a 
small M ej ). We take an example with L ~ 10 47 erg s^ 1 (B p ~ 10 14 G) and M ej ~ 1(T 4 M . 
To describe the dynamics in such a case, besides the spin down timescale T sd , we need three 
more characteristic time scales and the Lorentz factor value at the deceleration time 

T dec ~ 4.4 x 10 4 s ^XV" 1/Ifl 
T m ~ 3.6 x 10 3 s L^ 7 M cj „ 4 
T N2 ~ 4.5 x 10 7 s Ll%n- l ^TlH 

7 d ec ~ nZL^M^Jln-W + 1 (4) 

where Tm, Tm are the two time scales when the blastwave passes the non-relativistic to rela- 
tivistic transition line 7 — 1 = 1 during the acceleration and deceleration phases. With these 
parameters, one can characterize the dynamical evolution of the blastwave (FigJ2t). Based 
on the dynamics, we can quantify the temporal evolution of synchrotron radiation charac- 
teristic frequencies i/ a , u m , v c , and the peak flux, F V)mgx . The evolutions of the characteristic 
frequencies are presented in FigfSJb and collected in Table 1. 

At T dec , we have 

, a , dec ~ 5.0 x 10 8 Hz Ll^M^n^%\eZ, (^) (p + l) 3/5 /(p) 3/5 

z, m , dec ~ 1.3 x 10 14 Hz L%M^Lyi'\l_ x e l '% 
v c , dec ~ 9.6 x 10 14 Hz Ll%M^lln^e B z L\ 

F vmAec ~ 1.7 x 10 5 /iJy Ll'^Mf_ A n 2 ' 5 e l i%D^ (5) 

p/3p+22 Np/3p+2 n 

where f(p) = T , / 3p 1 + ig ^ / 3 pl u . With the temporal evolution power law indices of these param- 

>■( 12 - J 1 I 12 I 

eters (Table 1), one can calculate the X-ray, optical and radio afterglow lightcurves. Notice 
that there are two more temporal segments listed in Table 1, since u a crosses v m twice at 

T mal ~ 1.4 x 10 2 s L-^M![i^%%%]Ll ^ j (p + l) 1 ^) 174 , 
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T ma2 ~ 1.9 x 10 8 s Ll%n^T^eUe^ 2 (jz^)\p + l)^ f{p)~ 2/ \ (6) 

respectively. We present the light curves in X-ray (FigJ2E), optical and radio (10GHz) band 
(FigfSb). The distance is taken as 300 Mpc, the detection horizon of Advanced LIGO. 

Case II: M e j ~ M e j iC) i or T sd ~ T dec . The dynamics and the expressions of the charac- 
teristic parameters become simpler: 

T m ~ 12 s r'M e1 _ 4 T sdj3 



T N2 ~ 1.3 x 10 8 s e /3 M^lT Bd , 3 



cj,- 



7sd ~ 83.3eMT 1 _ 4 + 1 (7) 

The temporal indices of the evolutions of u a , u m , v c , F Utmax are listed in Table 1. 

As examples, we consider Lq ~ 10 49 erg s _1 (B p ~ 10 15 G) vs. M e j ~ 1O _4 M , which 
satisfies T sd ~ T dec . 

Similarly to Case I, we have 
, a , sd ~ 2.2 x 10 9 Hz e^M^/Xh^ (^) (P + l) 3/5 /(p) 3/5 

z, m ,sd ~ 2.7 x 10 17 Hz eM^_ A n^el_^_ 2 f j^|Y 
z/ c , sd ~ 8.6 x 10 14 Hz i-'M^_,n-^T- A %e B % 

F vmsA ~ 2.4 X 10 8 //Jy e 11 ^M9 M eT 8 -4^ 3/2 4 /2 -2^27 2 

/ _ 9 \-5/4 

T mal ~ i.4 x io- 1 s C^^n^^lUt* (fry) (p + i) 174 /^) 174 

T ma2 ~ 2.5 x 10 8 s ^^V%.<-A (j^f ) 2 (P + l)' 2/5 /(p)^ 2/5 (8) 

We also list the power- law indices for this situation in Table 1, and present the dynamics, 
typical frequency evolution, and the light curves in Figj3j 



Case III: M e j c l < M e j < M e j c 2 (T sd < T dec ). As an example, we take B p ~ 10 15 G, and 
ej ~ 10 3 M . 



For this example, the dynamics and the expressions of the characteristic parameters 
become 

T dec ~ 1.5 x 10 4 s £ V:i .W,VV' 13 



-8 - 



and 



T N1 ~ 59.9 s r^ej.-s^a 
T N2 ~ 2.7 x 10 7 s £ 1/3 n~ 1/3 

%d ~16.7£M-i_ 3 + l (9) 



, a , sd ~ i.6 x io 8 Hz (-h;;. v^.v), 5 , ( P + ip/bP 



1/2J2 .1/2 f P - 2 



z/ m , sd ~ 4.5 x 10 14 Hz rMTV 1/j! <-i6fl / ) -2 
z/ c , sd ~ 5.3 x 10 17 Hz C 4 Ml_ 3 n- 3 / 2 T- d 2 3 e B 3 !_ 2 2 



2 



JO — 1 

^, max ,sd ~ 6.5 x 10 2 /iJy ^M^n^T^e^D^ 2 

/ \ -5/4 

T mal ~ i.o s r 5/ X> 1/8l K^ e S5 (p + ^W 4 

T ma2 ~ 9.9 x 10 7 s e 1/5 n" 2/5 e 2 _ ie ^ /5 _ 2 ' (p + l)" 2 / 5 /(p)- 2/5 (10) 



The power-law indices of various parameters for this case are also collected in Table 1, 
and the dynamics, frequency evolutions, and light curves are presented in Fig. HI 

Case IV: M e j > M e j jCj 2. In this case, the blast wave never reaches a relativistic speed. 
The dynamics is similar to Case III, with the coasting regime in the non-relativistic phase. 
The dyn amics for a non - relativ istic ejecta and its radio afterglow emission have been dis- 



cussed in iNakar &: Piranl (120111 ) . Our Case IV resembles what is discussed in iNakar &: Piran 



(1201 lh . but the afterglow flux is much enhanced because of a larger total energy involved. 



3. Detectability and implications 

For all the cases, bright broadband EM afterglow emission signals are predicted. The 
light curves typically show a sharp rise around T sd , which coincides the ending time of the 



X-ray afterglow signal discussed by [Zhang ( 120121 ) due to internal dissipation of the magnetar 



wind. The X-ray afterglow luminosity predicted in our model is generally lower than that 
of the internal dissipation signal, but the optical and radio signals are much brighter. In 
some cases, the R-band magnitude can reach 11th at the 300 Mpc, if M e j is small enough 
(so that the blastwave has a high Lorentz factor) and the medium density is not too low. 
The duration of detectable optical emission ranges from 10 3 seconds to year time scale. The 
radio afterglow can reach the Jy level for an extended period of time, with peak reached 
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in the year time scale. These signals can be readily picked up by all-sky optical monitors, 
and radio surveys. The X-ray afterglow can be also picked up by large field-of-view imaging 
telescopes such as ISS-Lobster. 

Since these signals are originated from interaction between the magnetar wind and the 
ejecta in the equatorial directions, they are not supp osed to be a ccompanied with short 
GRBs, and some internal-dissipation X-ray afterglows ( jZhangl 120121 ) in the free wind zone. 
Due to a larger solid angle, the event rate for this geometry (orange observer in Fig.l) should 
be higher than the other two geometries (red and yellow observers in Fig.l). However, the 
brightness of the afterglow critically depends on the unknown parameters such as M e j, B p 
(and hence L ), and n. The event rate also crucially depends on the event rate of NS-NS 
mergers and the fraction of mergers that leave behind a massive magnetar rather than a 
black hole. 

This afterglow signal is much stronger than the after glow signal due to eje cta-medium 
interaction with a black hole as the post- merger product (INakar fc Piranll201l[ ). The main 
reason is the much larger energy budget involved in the magnetar case. Since the relativistic 
phase can be achieved, both X-ray and optical afterglows are detectable, which peak around 
the magn etar spindown time sc ale (10 3 — 10 5 s). The radio peak is later similar to the black 
hole case (INakar &: Piranll201ll ). but the radio afterglow flux is also much brighter (reaching 
Jy level) due to a much larger energy budget involved. The current event rate limit of > 350 
mJy radio transients in the minutes-to-days time scale at 1.4 GHz is < 6 x 10 -4 degree -2 yr -1 



(IBower fc Saull 1201 ll ). or < 20 yr all sky. In view of the large uncertainties in the NS- 
NS merger rate and the fraction of millisecond magnetar as the post- merger product, our 
prediction is entired consistent with this upper limit. Because of their brightness, these radio 
transients can be detected outside the Advanc ed LIGO horizon, w hich may account for some 
sub-mJy radio transients discovered by VLA (jBower et al.l 120071 ). 



Detecting these bright signals associated with GWB triggers would unambiguously con- 
firm the astrophysical origin of GWBs. Equally importantly, it would suggest that NS-NS 
mergers leave behind a hyper-massive neutron star, which gives an important constraint on 
the neutron star equation of state. With the GWB data, one can infer the information 
of the two NSs involved in the merger. Modeling afterglow emission can give useful con- 
straints on the ejected mass M e j and the properties of the post-merger compact objects. The 
combination of GWB and afterglow information would shed light into the detailed merger 
physics. 
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Table 1: Temporal scaling indices of various parameters in different temporal regimes for all 



dynamical cases. 
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Fig. 2. — Calculation results for Case I: L ~ 10 47 ergs -1 and M e j ~ 10~ 4 M Q (for all the 
examples, we adopt £ = 0.5, p = 2.3). (a) The dynamical evolution of the parameter 
(7 — 1); (b) Temporal evolutions of the characteristic frequencies z/ a , i/ m , and z/ c , and the 
peak flux F u ^ max ; (c) Analytical light curve in R-band (blue) and 10 GHz radio band (red); 
(d) Analytical light curve in X-ray band. In (c) and (d), the solid and dashed lines are for 
n = 1 and n = 10~ 3 , respectively. For n = 1 light curve, we mark the spectral and temporal 
indices for each segment. The main figures denote the time regimes when the light curves are 
detectable. The insets show the full light curves for completeness. Both X-ray and optical 
light curves reach their peaks around 10 4 s, and remain detectable in years. The radio light 
curve peaks around 10 7 s, and lasts even longer. The peak flux for X-ray, optical and radio 
could be as bright as 10~ n ergs _1 cm~ 2 , lOmJy and Jy, respectively. 
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Fig. 3. — Calculation results for Case II (L ~ 10 49 ergs -1 and M e j ~ 1O~ 4 M ). Captions 
are the same with Figure 1. For n — 1, both X-ray and optical light curve reach their peaks 
around 10 3 s, and the radio light curve peaks around 10 7 s. The peak flux of X-ray, optical, 
and radio is 10~ 9 ergs -1 cm~ 2 , lOOmJy and lOOmJy, respectively. Taking R-band magnitude 
20 and 10 -15 ergs _1 cm -2 as the detection limit, the durations of the detectable optical and 
X-ray afterglow are ~ 10 6 s and ~ 10 8 s respectively. The radio afterglow lasts even longer. 
For n = 10~ 3 , the signals for X-ray is still detectable, but with shorter durations, ~ 10 5 s. 



-14- 




Fig. 4. — Calculations results for Case III (L ~ 10 49 ergs -1 and M ej ~ 1O _3 M ). Captions 
are the same with Figure 1. For n — 1, both X-ray and optical light curve reach their peaks 
around 10 5 s, and the radio light curve peaks around 10 6 s. The peak flux for X-ray, optical 
and radio is 10~ 10 ergs -1 cm -2 , lOmJy and Jy, respectively. Taking R-band magnitude 20 
and 10 -15 ergs -1 cm~ 2 as the detection limit, the durations of the observable optical and 
X-ray afterglows are ~ 10 7 s and ~ 10 8 s respectively. The radio duration lasts even longer. 
If n = 10~ 3 , the optical signal just reaches R-band magnitude 20 around ~ 10 6 s, while the 
X-ray afterglow is detectable with a duration ~ 10 8 s. 
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